Existence of two-channel Kondo regime for tunneling impurities with resonant 

scattering 
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Dynamical tunneling systems have been proposed earlier to display a two-channel Kondo effect, 
the orbital index of the particle playing the role of a pseudospin in the equivalent Kondo problem, 
and the spin being a silent channel index. However, as shown by Aleiner et al. [Phys. Rev. Lett. 
86, 2629 (2001)], the predicted two-channel Kondo behavior can never be observed in the weak 
coupling regime, where the tunneling induced splitting of the levels of the tunneling system always 
dominates the physics. Here we show that the above scenario changes completely in the strong 
coupling regime, where - as a non-perturbative analysis reveals - the two-channel Kondo regime can 
easily be reached. We show that a tunneling systems end up quite naturally in this regime if the 
conduction electrons are scattered by resonant scattering off the tunneling impurity, and we also 
speculate about the possible origins of such a resonant scattering. 



I. INTRODUCTION 

There are a number of somewhat mysterious low- 
temperature transport anomalies in disordered point 
contactSji*2ii^ structurally disordered single crystals 
and disordered alloys ^^2^, which have not been ex- 
plained satisfactorily. In some cases the observed anoma- 
lies display power law behavioriS (~ T 1 / 2 or V 1 ^ 2 ), 
sometimes they show universal scaling properties^ but 
they have the common feature that all of them seem to 
be related to the presence of some dynamical impuri- 
ties. Although very attractive, electron-electron interac- 
tions fail to explain the absence of this zero bias anomaly 
in point contacts with strong static disorder^! and the 
complete absence of magnetic field dependence in some 
experiments^^ 

Much before the above-mentioned experimental re- 
sults, Vladar and Zawadowski proposed that dynamical 
two-level systems - abundant in amorphous regions^ - 
could lead to an orbital Kondo effect ^Zdl They consid- 
ered the motion of an ion in an effective double well po- 
tential, depicted in Fig. ^ interacting with free conduc- 
tion electrons. They assumed a simple potential scat- 
tering interaction between the atom and the conduction 
electrons, and derived the following effective Hamiltonian 





Fig. a 



Fig. b 



FIG. 1: Sketch of the tunneling system. The tunneling sys- 
tem is probably formed by a single atom in an amorphous 
region (Fig. a). The tunneling atom moves in an effective 
double well potential (Fig. b) . 



for low temperatures, where the atom moves by tunneling 
between the two minima of the potential well: 

„ A A z 

ti = t t T, 
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Here the Pauli matrices describe the motion of the 
particle in the double well potential, with r z = ±1 cor- 
responding to the left and right potential wells, Ao the 
tunneling, and A z the asymmetry of the potential (see 
Pig.HJ. The operators ip', o g create conduction electrons 
in some even/odd angular momentum channel (s and p, 
e.g.)r^ respectively, and are defined as 
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Here Do is a high energy cutoff discussed later and the 
^ {Oe/o o-'s satisfy canonical anti-commutation relations: 

This corresponds to a normalization of the fields ip e / 0iCr 
in Eq. JQ) such that the imaginary time propagators at 
T = temperature behave asymptotically as 
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The dimensionless couplings v x and v z in Eq. and 
the splitting A have also been estimated by Vladar and 
Zawadowski. Assuming a simple s-wave scattering U off 
the tunneling impurity they found 
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where go denotes the density of states of the conduction 
electrons at the Fermi energy, d is the tunneling distance, 
kp the Fermi momentum, Vb is the height of the tunnel 
barrier, and A the Gamow factor. The attempt frequency 
ujq is typically somewhat less than the Debye tempera- 
ture, and is typically in the range wo ~ 100 K assuming a 
tunneling atom of mass M ~ 50 m p , with m p the proton 
mass. Vladar and Zawadowski also obtained the pertur- 
bative scaling equations for the model defined by Eq. (JTJ 
and showed that the couplings and the dimensionlcss 
tunneling A = A /D satisfy the following differential 
equations 



dv x 
dl 



4VyV z 



dAo 
dl 



= (1- 



8v x {v 2 y +v 2 z ) , 
8(v 2 z +v 2 u ))A , 
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where I = \n(D /D) denotes the scaling variable with D 
the running cutoff (energy scale). The other equations 
are obtained by cyclic permutation from Eqs. © and J7J. 
The coupling v y is absent in the original Hamiltonian, it 
is, however, generated by the scaling procedure. Clearly, 
Eq. QBJ - apart from a prefactor in the second term - is 
that of the anisotropic Kondo model jii and generates a 
Kondo effect at the energy scale ~ Tk, the Kondo tem- 
perature, where the running coupling constants become 
of the order of unity. This Kondo effect is, however, strik- 
ingly different from the ordinary Kondo effect, where the 
low temperature physics is that of a Fermi liquid below 
Tk- In the present case spin is conserved in the course of 
scattering, and the two spin channels give rise to a singu- 
lar two-channel Kondo behavior below Tk, provided that 
the splitting of the two levels can be neglected^ The 
anomalous properties of this two-channel Kondo state 
(resistivity ~ T 1 / 2 , V/T scaling etc.) have been pro- 
posed to explain the experimentally observed anomalies. 1 
However, the tunneling amplitude Ao is always a relevant 
variable, and it ultimately kills orbital fluctuations (and 
thus the Kondo effect) of the tunneling system below a 
characteristic energy scale AJ. If this energy scale Aq is 
larger than Tk then the anomalous properties of the two- 
channel Kondo fixed point cannot be observed. Vladar 
and Zawadowski argued that - for a symmetrical tunnel- 
ing system - there is a parameter range where Aq <C Tk, 
implying that there is a temperature window where the 
physics is dominated by the Kondo effect. 

Unfortunately, Vladar and Zawadowski assumed in 
their analysis that the scaling equations above are valid 
at all energy scales below the Fermi energy Ep. How- 
ever, as pointed out by Aleiner et <d.~± this assumption 
is wrong: Electrons with energy £ 3> too follow the tun- 
neling particle adiabatically, 16 and do not give contribu- 
tion to the vertex renormalization. As a consequence, 
equations © and J7J) become trivial for too < D < Ep: 



for uo < D < E F . 



(8) 



dvi /dl = 
dh t /dl = 1 

Therefore, as pointed out by Aleiner et al., Eqs. |J?J| and 



(J7J) must be solved with the initial conditions D = ujq, 
A ~ e~ A > v x ~ (goU)(kpd) 2 e~ x . Since in the per- 
turbative regime Ao grows always faster that v x , Ag is 
always larger than Tk, and the two-channel Kondo ef- 
fect can never be observed for » z « 1, in contrast to the 
original conclusions of Vladar and Zawadowski. 

In this paper we first show that the arguments of 
Aleiner et al. hold only in the weak coupling regime, 
» z « 1. In the strong coupling regime, v z ~ 1 one must 
treat the coupling v z non-perturbatively and use the scal- 
ing equations originally derived by Vladar, Zimanyi, and 
Zawadowski in Ref. 20;, that treat v z exactly while han- 
dling the small coupling v x and the tunneling only in 
leading order. The analysis based on these equations 
clearly shows that there exist a critical value, v z = v ZyC = 
1/tt. Above this so-called Emery-Kivelson line^ v x typ- 
ically wins over Ao, i.e.,, the renormalized splitting is 
smaller then the Kondo temperature, Tk > Aq, and 
thus there is usually a wide temperature range where 
the physics of the tunneling system is dominated by the 
two- channel Kondo fixed point. 

Then we show that a possible and probably the most 
natural way to get into this regime, is to have a tun- 
neling atom that also acts as a resonant scatterer at the 
Fermi energy^ Indeed, as already observed in the orig- 
inal work of Vladar and Zawadowski, it is known from 
experiments that the value of v z can be large, suggestive 
of resonant scattering from the tunneling systems, what- 
ever they are. As we shall see, in this case it is not enough 
to start fro m the s imple potential scattering model stud- 
ied in Refs. Il2ll7l and one must take into account the 
dynamical motion of the internal levels of the tunneling 
atom. As we show later, due to this resonant scattering, 
one easily gets outside the range of validity of the per- 
turbative equations considered by Aleiner et al, Eqs. © 
and J7J), and ends up in the strong coupling regime. 

To show why resonant scattering is of primary impor- 
tance, let us consider a simple toy model describing an 
atom with a single resonant level at the Fermi energy, 



H 
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d a + h.c.) 



(9) 



Here d a is the annihilation operator describing the level, 
V is the hybridization, and creates a conduction elec- 
tron with spin a, momentum k, and energy £(k). This 
model can be trivially solved, and we can compute the 
spectral functions Qd(w) and q^uj) of the d-level and the 
fermionic field at the impurity site, ip a = X)k Ck >°'' 
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Here go is the density of states of the conduction electrons 
at the Fermi energy, and T = 2irg V 2 is the width of the 
resonance. Clearly, if we now try to move the atom, it will 
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FIG. 2: Density of states (spectral function) of a resonant 
level d and the conduction electron operator ip hybridizing 
with it. The conduction electrons' spectral function is sup- 
pressed While Qd displays a huge resonance at the Fermi en- 
ergy. 



then couple to its own resonant level d the most strongly, 
which is centered at the atom, and not to the rest of 
the conduction electrons, which are sitting at neighbor- 
ing ions in reality. Even for a not too narrow resonance 
with r ~ 1 eV, the density of states of the resonant level 
at the Fermi energy is large, ~ 1/T, while that of the con- 
duction electrons hybridizing with the atom actually van- 
ishes at the Fermi energy (see Fig. - This large increase 
in the density of states is what immediately pushes the 
system in the strong coupling regime, v z > v z%Cl where 
the two-channel Kondo behavior prevails, provided that 
the resonance is narrow enough. For typical parameters 
we obtain that, the resonance should be narrower than 
about 0.1 E F - 0.01 E F ~ 1000 - 10000 K. 

We emphasize though that the above resonance does 
not have to be narrower than luq: As we discussed earlier, 
only electrons with energy \uj\ < luq are unable to follow 
the motion of the atom at a time scale ~ l/u>o- These lat- 
ter are exactly the electrons that are responsible for the 
Kondo effect. However, these electrons only need a time 
scale ~ 1/r to notice the increased scattering strength 
at the Fermi energy, i.e., all these 'slow' electrons see an 
increased scattering strength from the atom. This simple 
picture can be readily verified^ through a path integral 
treatment similar to that of Ref. as well as by the ex- 
tension of the the diagrammatic approach of Refs. IT^ol 
For a resonance narrower than luq our analysis must be 
slightly modified, and the vertex renormalization only 
occurs in the range D < T, since electrons with energy 
lu > r do not couple strongly to the atom. 

As we show later, the increase in the coupling constant 
can simply be understood as a matrix element effect: the 
wave functions of the conduction electrons at the Fermi 
energy (where the resonance is located) have an increased 
amplitude at the tunneling atom's position, and therefore 
the motion of the atom couples more strongly to the con- 
duction electrons. 



It is thus the structure of a resonant level itself which 
is ultimately responsible for the orbital Kondo effect in 
the above scenario, and this resonance cannot be re- 
placed by a simple potential scatterer in a free electron 
gasi^i The technical reason for this is that in a sim- 
ple potential scattering model there is always a term 
~ 2 ff ^e.o^o./^jta^uo-! which is large, cannot be ne- 
glected, leads to a suppression of the conduction elec- 
tron's density of states, and ultimately reduces the effec- 
tive value of v z . As we shall see, the effective Hamiltonian 
for this resonant level finally takes also the form, Eq. QJ, 
however, the couplings will be related to some atomic 
orbitals, and no large term ~ vo u, ' t l ) }ia' l l } v a PP ears - 
Therefore we can use the non-perturbative scaling equa- 
tions valid for all values of v z , derived longtime ago by 
Zimanyi, Vladar, and Zawadowski in Ref.f20|to construct 
the phase diagram of our model. The summary of this 
analysis is shown in Fig. OH The bell-shaped line in Fig. 
shows the estimated value of the Kondo temperature, 
Tk, while Ag is the renormalized tunneling amplitude. 
In Fig. we also show another energy scale, A|, associ- 
ated with a two-electron scattering process^ Below this 
energy scale the two-channel Kondo behavior is also sup- 
pressed. As one can see from Fig. 01 there is an extended 
regime dominated by the two-channel Kondo fixed point. 
The size of this region increases if we decrease the value 
of v x and Ao, however, then Tk is also shifted towards 
smaller temperatures. As we discussed above, the tunnel- 
ing system ends up immediately in this regime if there is 
a sufficiently narrow resonant scattering on the tunneling 
atom. 

We would also like to mention that another possibility 
to get to the strong coupling regime is to increase the 
tunneling distance and have a broader resonance at the 
Fermi energy. This might be possible for very light tun- 
neling impurities, possibly by hydrogen: hydrogen has 
typically resonant scattering in the s-channel and it can 
possibly tunnel over a large distance in a metal. Zero bias 
anomalies have indeed been observed in hydrogen doped 
palladium point contacts where the anomaly is clearly 
related to the presence of hydrogen^ 

The rest of the paper is structured as follows: In Sec- 
tion HU we shall construct a microscopic model that en- 
ables us to treat the resonance properly and map this 
problem to the original Vladar-Zawadowski model. We 
then analyze the scaling equations and construct the 
phase diagram of this model in section IIIII Finally, in 
the concluding section we shall speculate about possible 
candidates producing a resonant scattering at the Fermi 
energy. 



II. MICROSCOPIC MODEL 

As a first step in constructing a microscopic model for 
a tunneling system coupled to the electrons, we write the 
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FIG. 3: Phase diagram of the tunneling system as a func- 
tion of the phase shift S = artan(u z 7r). The arrow indicates 
the two-channel Kondo regime. The symbol Aq denotes the 
energy scale where the tunneling becomes dominant, while 
A2 is the energy scale where a special two-electron scatter- 
ing process becomes important. We have used u>o ~ 100K, 
Do ~ 10 5 K, Ao = 10K, and v x « 0.1 as bare parameters. 



Hamiltonian in a first quantized form as 
1 



H 



Afi 

2M ^ 2m 



Y—A: 



(12) 



c/(r) + e n*i -r)+EE v ^ - R 

t j#0 i 



■31 ' 



where R denotes the coordinate of the tunneling atom, 
and rVs are the coordinates of the conduction electrons. 
We treat the motion of the tunneling atom quantum me- 
chanically, however, we assume that all other atoms are 
immobile, and their position Rj (j = 1, . . .) is a con- 
stant classical variable. The tunneling atom has a mass 
M and moves in the double well potential C/(R) formed 
by the rest of the ions. In the Hamiltonian H12|) we also 
assumed that electrons form a Fermi liquid and thus ne- 
glected electron-electron interaction, which is supposed 
to be included in the effective electron-ion interaction 
potential V(r — R) at the Hartree or Hartree-Fock level. 

To simplify this Hamiltonian, as a first step, we rewrite 
Eq. O as H = H + ff lnt with 



Hn = 



i j—0.... i 

E (V(Ti - R) - V(r. 



R-i) , 
R )) , 



(13) 



(14) 



justified by the observation that the tunneling impurity 
moves in a small cavity, our analysis does not rely on it, 
and our conclusions are independent of this approxima- 
tion. To arrive at a tight-binding Hamiltonian, we first 
solve the atomic Schrodinger equation: 



-2^ + F(r - 



R o) M r ) = e P M r ) 



(15) 



where labels the atomic levels. The states ip^ above 
are atomic eigenstates centered at the tunneling atom. 
As a next step, we solve the Schrodinger equation for 
the rest of the conduction electrons without making any 
approximation, 



f?n(r) = U Vn{r) ■ (16) 



Then, using the wave functions ip n and r/„, we can com- 
pute the appropriate overlap matrix elements and con- 
struct the following tight binding Hamiltonian for the 
conduction electrons: 



+ E (W^A* + h - 



(17) 



where c n ^ denotes the creation operator of a conduction 
electron in the (extended) state n with spin a, a cre- 
ates a conduction electron at the atomic orbital fi, and 
tn,p, denotes the corresponding hopping matrix element. 
Since the interaction part Hi n t is large only at the posi- 
tion of the tunneling impurity, we can integrate out all 
electrons c n a and arrive at the following imaginary time 
effective action for the d-levels: 



D eff — 



E f dr f dr'd^g-i (r - r')^v(r') 



(18) 

where Q ufl i denotes the local propagator of the ci-levels, 
which can be expressed in Fourier space as 



i 



n,u • /■ "7l,/r 

' p lOJ n - t,n 



(19) 



This Green's function can be written even more conve- 
niently in a spectral representation as 



where Ro is a somewhat arbitrary 'reference position' 
of the ion, which is chosen to minimize the interaction 
part H mt . A natural choice is, of course to choose Ro 
to correspond to the maximum of the barrier in J7(R). 
Clearly, the non-interacting part H can be diagonalized. 

To make further progress, we shall adopt a tight- 
binding scheme for the atom. Though this approach is 



(20) 



It is this spectral function Qua' (^) which contains the 
resonance discussed in the previous section, and which 
solely determines the low temperature properties of the 
tunneling system. 
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In the following we shall consider the simplest case, 
where the tunneling system is fully symmetrical. We first 
diagonalize the Hamiltonian of the tunneling atom, 



(21) 



In principle, we could formulate our theory by keeping 
all levels of the tunneling particlej2i However, at low 
enough temperatures only the two lowest lying even and 
odd states, $ e and <i> matter, and the role of all other 
eigenstates is to reduce the electronic cutoff to a value of 
the order of the Debye frequency, Do — > ui^M^ There- 
fore, in the following we shall keep only these two states. 
To obtain the usual tunneling form that occurs in Eq. Q , 
we introduce the left and right states & T =± 



1 



V2 



(22) 



which transform into each-other under reflection, and 
rewrite the Hamiltonian Eq. I|21|l in this restricted sub- 
space as 



rrtun _ A 

H a - — 2 x ' 



(23) 



with t x the Pauli matrix. The tunneling amplitude A = 
E Q — E e is approximately given by Eq. @. 12 

Having diagonalized the non-interacting part i|13|) of 
the Hamiltonian, we now turn to the analysis of the in- 
teraction part. First, we shall simplify our treatment by 
keeping only those two even and odd electronic states 
(p = e,o), which couple the most strongly to the tunnel- 
ing particle. In a real system, these correspond to the 
most strongly coupled d and p-states, or s and p states. 
In case of a p-state, e.g., it is clear that the p-state aligned 
along the tunneling axis couples the most strongly to the 
motion of the atom. Due to the assumed spatial symme- 
try of the Hamiltonian, the corresponding spectral func- 
tion is diagonal in the index p, q^i = S^^g^. We shall 
assume furthermore that at least one of the spectral den- 
sities, say g contains a resonance of width T, while the 
other one is approximately constant, g e sa go ~ 1/Ep 



9o{u) = 

Qe{u) = 



1 



2tt lu 2 
Qo ■ 



T 2 /4 



(24) 
(25) 



To construct the interaction part of the Hamiltonian in 
this restricted basis, we have to compute matrix elements 
of fli n t given by Eq. (|14(l . Assuming that the tunneling 
distance is small, we can approximate fli n t as 



fli 



q=x,y,z 



R, 



dV 



q 3r 



- Y R p R q 



d 2 V 
dr p dr q 



+ 



(26) 

and then compute the appropriate matrix elements to 
obtain the effective Hamiltonian as 

Hut « A(|$ e )(* | + \®o){$e\)(di a d oa + dt a d e(T ) 
+ (B e \$ e )($ e \ + B \<P )($ \)(C e dld ea + Codi^) , 



where summation is assumed over repeated spin indices, 
and the constants A, B e / a and C e / are given by the 
following integrals: 



.4 



B. 



/o 2 \ e/o 



dr z 

$ K /n\Rl\$ 



d 2 V(r) l 



e/o! ) 



We/o) 



(27) 
(28) 
(29) 



This Hamiltonian can easily brought to the form, Eq. QJ, 
if we notice that only electronic excitations with ener- 
gies < loq contribute to the orbital Kondo correlations. 17 
Therefore, if the width of the resonance is broader than 
the attempt frequency, r > ojq ~ WOK ~ O.OleV, then 
we can introduce the new fermionic fields normalized ac- 
cording to Eq. (J2J) 



dea/ \fgl i d oa ~ > Ipoa = d oa / ^fg~ , (30) 



where g e and g denote the density of states at the Fermi 
energy in the even and odd channels, respectively. In 
terms of these fields, the interaction Hamiltonian takes 
on the form Eq. ft}, with the couplings v x and v z given 
by: 



v z = A^g e g Q , 

v x = j(fl e - B )(C e g e - C g ) 



(31) 



There is two more terms that appear in addition to these 
two terms, both of which are small: One is a simple 
potential scattering term, ~ (V4V"e — iplrfo) that only 
renormalizes the spectral densities g and g e , and can 
be eliminated by a simple counterterm procedure. The 
other term is proportional to ~ T x (^pl^p e + iplipo), and 
gives a small renormalization of the double well poten- 
tial U (R) This term can be taken into account at the 
Hartree level and is of no importance. In the following 
we shall therefore keep only the two couplings v x and v z 
in Eq. (ED- 

Note that the above procedure of integrating out the 
electrons c„ i(T and then rescaling the couplings is just a 
technical trick to extract the dimensionless couplings and 
to derive quickly the scaling equations. However, one can 
proceed in the usual way as and obtain the dimensionless 
couplings by analysis the str ucture of the dimensionless 
vertex functions, as in Refs. |25l33| . 

The couplings v 7 , and v x can be estimated along similar 
lines as in Refs. [i2lT5t and one obtains: 



y/g e QoWe\d— \<p ) 



-A — — g (>f \d 2 



A_J_ 

'16 Vb' 



{dr z f 



<Po) 



(32) 
(33) 



where in the second equation we assumed a quartic dou- 
ble well potential with barrier height Vb, and displayed 
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only the contribution of the odd channel. The constant 
d in the previous equations denotes the tunneling dis- 
tance. A similar but smaller contribution is given by the 
non-resonant even channel. Remarkably, the matrix ele- 
ments above are expressed in terms of the atomic (tight 
binding) orbitals of the tunneling atom. They can be 
easily estimated assuming a simple Coulomb interaction, 
V(r) = — e 2 /r, and hydrogen-like wave functions. Taking 
the Is and 2p orbitals, e.g., we find 



(2p 



dr z 
d 2 V{r) 



2p) 



ao 21\/2 a 
e 2 1 



a 40 ag 



5.4 



-2.43 



eV 

T ' 

eV 
12 



(34) 
(35) 



where ao is the Bohr radius. 

It is not difficult to show that even for a wide resonance 
with T - 1000 - 10000 K and d ~ 0.3l, v x is about the 
same as the dimensionless tunneling Ao /ujq. More impor- 
tantly, however, v z is proportional to l/\/r. Therefore, 
if the conduction electrons scatter resonantly off the tun- 
neling atom, then the value of v z can become large, and 
eventually become larger than the critical value corre- 
sponding to the Emery-Kivelson line mentioned in the 
introduction, u 2jC = l/it. For typical parameters this 
transition takes place where the width of the resonance 
is around T ~ 1000 - 10000 K. 

Throughout this paper, we shall assume that r > luq, 
i.e., that the resonance is relatively broad compared to 
typical frequencies related to atomic motion. Under these 
conditions, the cutoff energy is given by loq and the cou- 
plings between the tunneling system and the electrons 
can be safely approximated by their values at the Fermi 
energy. 

The above large increase in the effective couplings can 
be understood as a simple matrix element effect, and is 
related to the well-known structure of scattering states: 
To clarify this point, let us consider the textbook example 
of scattering states in the s-channel of a simple spheri- 
cally symmetric potential scatterer, U(r) = V 5(r — ro). 
The scattering wave functions in this simple case take 
the form sin(fcr + S)/r and b sm(kr)/r for r > r$ and 
r < ro , respectively, with k the radial momentum of the 
electrons and 6 the scattering phase shift. The amplitude 
b 2 of the wave function inside the sphere is given by the 
expression, 



b 2 = 



k 2 



k 2 + V 2 sin 2 (yfc) + 2kV sin(2fc) ' 



(36) 



Where we used units of ro = 2m = h = 1 . For V = the 
amplitude of the wave function is simply one. However, 
for larger values of V a resonance appears, and b 2 displays 
a sharp peak as a function of energy (Fig.EJ. It is a trivial 
matter to show that the amplitude b of the resonance is 
simply related to the width T of the resonance, 6 ~ 1/ y/T, 
and becomes larger and larger for sharper and sharper 
resonances. This large factor b shows up in the local 




FIG. 4: Amplitude b of the scattering state of the conduction 
electons at the origin in a simple scattering state model. The 
amplitude of the scattering state at the origin increases as 
the resonance gets sharper and sharper. The inset shows the 
potential producing the resonance. 



density of states and also any matrix element computed 
in terms of the appropriately normalized scattering states 
and results in an increase of all couplings, provided that 
the resonance appears at the Fermi energy. 



III. NON-PERTURBATIVE SCALING 
ANALYSIS 



We have seen in the previous section that for resonant 
scattering from the tunneling impurity the coupling v z 
can be very large. In this case the perturbative scaling 
equations, Eqs. (I(il7ll are insufficient and a new approach 
is needed. Fortunately, for Ao, v% -C 1 one can construct 
scaling equations which are non-perturbative in the cou- 
pling v z by generalizing the computations of Yuval and 
Anderson. 20-22 In this limit one obtains the following scal- 
ing equations 



c ,20 



oA 

dl 

dv x 
~dl 

dv. 



1 — 81 — 

7T 



Ao 



.5 

i—v. 



TV 



y 



(37) 
(38) 
(39) 



with the phase shift S defined as 5 = arctan(7rw z ) . Here 
we neglected terms of order 0(v 2 , v 2 , Ag). These terms 
give rise to a renormalization of 5 in the strong coupling 
regime, and slightly change the numerical values of the 
various energy scales we estimate, but do not affect their 
overall scale and the picture obtained using the above 
equations. Equations l|38|l and (|39|l can be rewritten in 



terms of the average coupling, v± 



y )/2, and the 
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asymmetry, V- = (v x — v y )/2 as 



dv± 
~dT 

dv- 

~dT 



7T 



4* 

7T 



(40) 
(41) 



In other words, the anisotropy v x ^ v y is irrelevant, while 
the average coupling v± is relevant. Therefore the two 
couplings v x and v y become rapidly equal in the initial 
stage of the scaling. Eqs. (|40|l and l|41|) are only valid 
below the scale D w ujq; above this energy scale the 
couplings v x , v y remain unrenormalizcd, while Ao trans- 
forms according to its engineering dimension and satisfies 
Eq. JHJ). Therefore Eqs. (flTij) and l|IT]l must be solved with 
the initial conditions D — ojq, Ao(^o) = Aq/ojq ~ e~ A , 
v y (u>o) = and v x (ojq) given by Eq. 

From Eqs. lET7|) and we see that the value v z = l/n 
corresponding to 6 — 7r/4 is very special: For S < 7r/4 
the scaling dimension of Ao is always larger than that 
of v x , and therefore tunneling is always more relevant 
than assisted tunneling. As a result, the motion of the 
tunneling system freezes out at an energy scale Ag <C Aq 



AS = AS,> = A„ ( ^° 

1 Wo 



8 



i-8 (ir 



(42) 



where the assisted tunneling is still small compared to 
unity, and thus the system is far from the two-channel 
Kondo fixed point. 

Eq. I|42|) only makes sense if Ag > Tk- This is always 
the case for 6 < ir/4. If, however, S > tt/A then the 
physics can become very different. In this case v x grows 
faster than Ag, so v x ~ 1 can be satisfied first. The con- 
dition v x ~ 1 defines the so-called Kondo temperature: 



TV 



(43) 



where the renormalized tunneling satisfies 



A (Tj 



K 



2A 



K~ 



UJ 



JL ( v -± 

v x \2 



4» _ i 

7T 

4A (i _ 2i) 

7T V 7T / 



(44) 

From this equation it is obvious that for 5 > 7r/4 the 
effective tunneling amplitude is usually still small at the 
energy scale T ~ Tr- compared to the temperature itself. 
Below the Kondo energy Tk the scaling equations are 
governed by the two-channel Kondo fixed point. There 
the tunneling is still relevant and has a scaling dimension 
1/2 and therefore its scaling equation must be replaced 
byis 



dA 



(D < T K ) 



(45) 



Integrating this equation we obtain for the renormalized 
tunneling amplitude for the case Ag < Tk, 



A 



o,< 



T K 



( 2A ' 



2ii 



<^T K 



(46) 



From our discussions immediately follows that in the 
tunneling regime for a symmetrical tunneling center with 
v z > v z , c there is typically a large temperature win- 
dow, Aq < T < Tk, where the two-channel Kondo 
fixed point rules and the non- Fermi liquid properties such 
as the ~ s/T resistivity anomaly associated with the 
two-channel Kondo fixed point should be manifest. We 
have to emphasize, however, that below the scale Ag the 
physics becomes again that of a boring Fermi liquid. The 
corresponding crossover lines were sketched in Fig. [31 for 
typical parameters of the tunneling system. As one can 
see in the figure, there is a large region in the parame- 
ter space, which is governed by the two-channel Kondo 
behavior. This region increases even further for smaller 
values of Ao ~ v x , however, it also shifts towards smaller 
values of T K - 

Before we conclude this section, let us discuss an- 
other important issue, raised by Moustakas and Fisher. 23 
Moustakas and Fisher observed that at the two-channel 
Kondo fixed point a special two-electron scattering of the 
form [V>± = (^ e ±Vo)/V2], 



H 2 = -^r x (V'+ T V , + i V'-|V'-T + h - c ) 



(47) 



is also relevant, and that for 8 > 7r/4 this operator is more 
relevant than the tunneling Ao discussed above. The 
coupling A2 above denotes a dimensionless coupling con- 
stant, and its bare value can be estimated to be around 
v x , while the energy scale Dq is of the order of the Fermi 
energy Ep. Below the energy scale D = uj this operator 
satisfies the scaling equation: 

2 



dA 2 
dl 



l + 4~ 8 (^) ) A,. (4N, 



Fortunately, this operator is irrelevant at high energies, 
u>o < D, where its scaling dimension is simply —1, 



dA 2 
dl 



(D > u ) . 



(49) 



As a result, by the time we get into the regime D < ujq 
this two-electron process is reduced by a factor u>q / Dq ~ 
1CP 3 compared to v x . 

Below ujq this process becomes relevant, and generates 
a cross-over to a Fermi liquid at an energy scale Aj . Sim- 
ilar to Ag, we have to distinguish two possibilities. For 
larger values of 6, Tk < A|, and A 2 is given by the 
following expression: 



2,> 



w 



^A 2 
E F 



= ^- 8 (f, -1 

7T \ 7T . 



(50) 



(51) 
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However, this formula is not correct if Tk > AJj, 
since A 2 has scaling dimension 1/2 below the Kondo 
temperature^ Carrying out an analysis similar to the 
case of Aq < Tk we find in this regime that 




The corresponding crossover line is also shown in Fig. [31 
As one can see, A 2 does not play a significant role in the 
regime where Tk is the largest, however, for larger values 
of 5 it is indeed A 2 that provides an infra-red cutoff for 
the two-channel Kondo behavior. 

IV. CONCLUSION 

In the present paper we constructed a theory which 
describes the low temperature behavior of a tunneling 
particle with resonant scattering. We have shown that 
the resonance plays a crucial role in the physics of the 
tunneling center: While the local density of states of the 
resonant atomic state has a huge peak inversely propor- 
tional to the width of the resonance, Qd ~ 1/r, the den- 
sity of states of the corresponding conduction electrons 
is strongly suppressed, and thus the tunneling atom cou- 
ples the most efficiently to its own electronic state. De- 
spite this complication, after a proper treatment of the 
electron-tunneling particle interaction, we arrive at the 
original model of Vladar and Zawadowski. However, our 
treatment does not rely on the free electron approxima- 
tion, the various coupling constants can be large and they 
are determined by atomic integrals. 

As we have shown, the resonant scattering can com- 
pletely change the physics of the tunneling impurity, and 
push it in the non-perturbative regime. While in the 
weak coupling regime the tunneling always intervenes 
before we enter the two-channel Kondo regime, in this 
non-perturbative regime the two-channel Kondo behavior 
usually dominates over the spontaneous tunneling over 
a wide temperature range of possibly several decades for 
physically relevant parameters. We emphasize again that 
a model where the tunneling center interacts with some 
local interaction with a free electron gas is unable to cap- 
ture this behavior, iL24 The physical reason for this is 
simply that electrons are not free. Even if we do not 
move the tunneling atom, the electron's wave function is 
adjusted to the atomic potential of the tunneling impu- 
rity. This effect has been heuristically taken into account 
in Ref. I2M where only the change in the scattering po- 
tential has been considered as a perturbation, without 
justifying this approach. Also, the simple potential scat- 
tering model is unable to take into account the dynamics 
(retardation effects) of the electronic states at the tunnel- 
ing impurity itself, which one actually tries to eliminate 
from the theory. 



It is not quite clear what the origin of such a resonant 
level could be in practice. In the point contact measure- 
ments, one natural candidate would be hydrogen. Hydro- 
gen is small enough, can diffuse into the substrate, and is 
known to have a phase shift 7r/2 in the s-channel* 6 [This 
phase shift, which just characterizes the atomic scatter- 
ing off a Hydrogen ion is not to be confused with the 
phase shift S corresponding to v x .] Moreover, hydrogen 
is light: This implies that the cutoff loq can be in the 
range of ujo ~ 1000 K. As a result, all temperature scales 
can be about an order of magnitude larger for hydro- 
gen than the ones in Fig. |21 and Tk can be easily in 
the range of a few Kelvins. Being light, hydrogen can 
tunnel over relatively large distances, d ~ 1 1, imply- 
ing that v z is presumably large even if the resonance is 
very broad. In fact, in the point contact experiments it is 
very difficult to exclude the presence of hydrogen in the 
course of sample preparation^ and recent experiments 
on hydrogen doped palladium point contacts indeed ex- 
hibit zero bias anomalies associated with the presence of 
hydrogen. 28 However, more detailed calculations would 
be needed to estimate the size of the coupling v z for tun- 
neling hydrogen. 

Tunneling systems with small effective masses can be 
formed by dislocations tooiSS In this case, however, the 
spatial extent of the defect can be large, and it is not quite 
clear how the two orbital scattering channels driving the 
two-channel Kondo effect could be selected. 

Other natural candidates would be impurities with 
strong magnetic correlations, i.e. Kondo- like impurities. 
In this scenario, two types of Kondo effect take place: 
a magnetic Kondo effect with a large Kondo tempera- 
ture, T^ agn , and the orbital Kondo effect discussed so 
far at a much smaller energy scale, Tg- rb <C T^ agn . The 
magnetic Kondo resonance would provide the resonance 
needed for the orbital Kondo effect, and the magnetic 
correlations serve only to boost up the couplings of the 
tunneling system and generate an orbital Kondo effect 
at lower temperatures. Above the magnetic Kondo tem- 
perature, T > Tj£ asn , magnetic scattering provides also 
a strong inelastic scattering. Correspondingly, there is 
a temperature-dependent time scale, r sp i n (T) at which 
the spin of a conduction electron is flipped. These spin- 
flip processes could, in principle, destroy the two-channel 
Kondo behavior. To be able to neglect these spin flip pro- 
cesses one needs to satisfy the criterion l/rfli p (T) < T at 
all temperatures. Luckily enough, the rate 1/tq\ v {T) is 
suppressed below the Kondo scale. Unfortunately, how- 
ever, to our knowledge, this spin flip rate has never been 
determined so far reliably. However, one can obtain a 
simple estimate for it using the knowledge about the in- 
elastic scattering cross section^ From these considera- 
tions one concludes that well below Tk, in the Fermi 
liquid range the spin flip rate must scale as ~ T 2 /Tk, 
while in the vicinity of Tk it must be of the order of 
lA'spin(T') ~ T. Above Tk this rate must scale to zero 
logarithmically. From these considerations it appears 
that the criterion l/Vfli p (T) < T is always satisfied, and 
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one can therefore probably safely neglect spin flip scatter- 
ing processes. We have to emphasize though that these 
impurities do not need to be Kondo impurities in the 
usual sense, since the width of the magnetic 'Kondo res- 
onance' can be in the range of thousands of Kelvins or 
even larger. For such a correlated impurity the local den- 
sity of states remains unrenormalized. 30 However, there 
is a strong field renormalization proportional to the Z- 
factor, Z ~ T™ as " /Ep, which ultimately rescales v z as 
v z —> Vg/\Z, We emphasize again that even an ex- 
tremely large (magnetic) Kondo temperature in the range 
of ~ 1000— 10000K could give rise to the phenomena dis- 
cussed in this paper, and therefore the usual transport 
and specific heat anomalies associated with the magnetic 
Kondo effect may be hardly observed in this case. 

In the case of a tunneling system with strong mag- 
netic correlations, the zero bias anomaly may be very 
sensitive to the external field too. The reason is that in 
this case the density of states and thus the couplings in 
the two spin channels may depend rather sensitively on 
the applied magnetic field. This translates to a chan- 
nel anisotropy in the effective two-channel Kondo model 
and drives the system to a Fermi liquid. This may possi- 
bly explain the strong magnetic field dependence in some 
experiments^ 

Strongly coupled electron-phonon systems have also 
been proposed as possible candidates to produce an 
orbital Kondo effect^ In this case the two- level sys- 
tems form dynamically. It is, however, not clear how 
the regime of extremely strong electron-phonon coupling 
needed can be reached. Furthermore, in these studies 
only a few vibrational modes could be considered, which 



has been shown to be insufficient to produce the correct 
low-temperature behavior ^fiiSS i 

Finally, let us comment on the presence of the splitting 
A z . In our previous analysis we completely neglected the 
asymmetry of th e tun neling centers. In lattice structures 
such as in Refs. |7l9lll0l ] A z can be rather small. How- 
ever, in a disordered point contact A z has a random dis- 
tribution and should provide a low energy cutoff for the 
non-Fermi liquid properties similar to the spontaneous 
tunneling. Therefore a break-down of universal scaling 
is expected due to the presence of A z . Such a break- 
down of universal scaling has indeed been observed in 
Ti point contacts^ where the zero bias anomaly did not 
depend on the presence of an external magnetic field, 
was sensitive to electro-migration, and had an amplitude 
consistent with the presence of just a few tunneling cen- 
ters. These experiments seem to be in perfect agreement 
with all predictions of the two-level system model. How- 
ever, it is much harder to explain the origin of the zero 
bias anomaly in Cu samples »i While the sensitivity to the 
magnetic field could be explained assuming that the tun- 
neling impurities have a sharp resonance, the zero bias 
anomaly in these experiments has a very large amplitude, 
and as pointed out by Smolyarenko and Wingreen;2i the 
observed anomalous resistivity exponent close to 1 /2 can 
hardly be understood assuming a completely random dis- 
tribution of A z . 

I am grateful to B. Altshuler, I. Aleiner, L. Borda, Sz. 
Csonka, A. Halbritter, and A. Zawadowski for valuable 
discussions. This research has been supported by Hun- 
garian Grants Nos. T046267, T046303, and T038162. 
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I. e., the atom itself that moves has an electron structure 
with resonant scattering at the Fermi energy, even if its 
position is fixed. 



